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Abstract 

The contribution of the quark self energy to the meson masses is studied in the framework of the 
relativistic flux tube model. The equivalence between this phenomenological model and the more 

o 

QCD based rotating string Hamiltonian is used as a guide to perform the calculations. It is shown 
that the addition of the quark self energy to the relativistic flux tube model preserves the linearity 



Q 

co 



PACS numbers: 12.39.Pn, 12.39.Ki, 14.40.-n 



of the Regge trajectories. But, following the definition taken for the constituent quark masses, the 

Regge slope is preserved or decreased. In this last case, experimental data can only by reproduced 

by using a string tension around 0.245 GeV 2 . Two procedures are also studied to treat the pure 

\Q . flux tube contribution as a perturbation of a spinless Salpeter Hamiltonian. 
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I. INTRODUCTION 



A successful way of understanding the properties of the mesons is to approximate the 
gluon exchanges between the quark and the antiquark by a string (the QCD string), char- 
acterized by its energy density, or tension. The relativistic flux tube model (RFTM) is a 
phenomenological model based on this picture 

HQ. 

More recently, the rotating string 
model (RSM) has been derived from the Nambu-Goto Lagrangian as an effective model 

nn 

which also describes a meson as a quark and an antiquark linked by a string |3, |4J. The 
physical content of these models is very similar, and it has been shown that they are actually 
classicallv equivalent once the auxiliary fields appearing in the RSM are completely elimi- 
nated 0,0]. A rather good description of the experimental meson spectrum can be obtained 
with the original RFTM supplemented by a Coulomb term []]]. But unfortunately, a strong 
coupling constant larger than what it is expected from experimental analysis and lattice 
simulations must be considered. The necessity of finding new contributions arising from 
neglected physical mechanisms is thus clear. Recently, a quark self energy (QSE) contribu- 
tion was introduced, which is due to the color magnetic moment of the quark propagating 
through the vacuum [^J. The QSE brings a negative contribution to the hadron masses, 
and seems to be an interesting way of reproducing the experimental data with a smaller 
Coulomb term. 

Our purpose here is to study the influence of the QSE on the meson masses, especially 
using the RFTM. Our paper is organized as follows. The Sec. |TI] is a short presentation of 
both classical RFTM and RSM, where we also underline their classical equivalence. Since 
the equations describing these models are quite complicated, approximate equations are 
developed in Sec. IHIl from which an analytical mass formula is derived in Sec. IIVI In 
Sec. El the quality of the approximate equations are compared with the original quantized 
RFTM. We then introduce the QSE in Sec. I VII and discuss its effects on the meson spectrum 
in Sec. IVHI A comparison with experimental data is performed in Sec. IVIIII Finally, some 
concluding remarks are outlined in Sec. IIX1 
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II. ROTATING STRING AND RELATIVISTIC FLUX TUBE 



We will present in this section two effective meson models : the rotating string model 
(RSM) and the relativistic flux tube model (RFTM). Our purpose is to give a presentation 
of the principal features of these models, and to underline their classical equivalence. 

It has been shown in Ref. [3[ that starting from the QCD Lagrangian, the Lagrange 
function of a meson with spinless quarks can be built from the Nambu-Goto action. For two 
quarks with masses mi and m,2, and a string with tension a, this action has the well-known 
form 

C = — m\\J~x\ — m 2 y^ — a J d(3\J (ww') 2 — w 2 w' 2 . (1) 

In this action, Xi is the coordinate of quark i and w is the coordinate of the string, w depends 
on two variables defined on the string worldsheet : one is spacelike (3 and the other timelike 
t. We have also defined w' = dpw and w = d T w. Introducing auxiliary fields (also known 
as einbein fields) to get rid of the square root in and making the straight line ansatz to 
describe the string connecting the quark and the antiquark, an effective Hamiltonian can be 
derived, which reads 

1 \p 2 + m 2 vl + m 2 2 

H(fli, V) = 1 h fl X + fjQ 

2 L A*l 

+a 2 r 2 — + df3u ■ 1 



(2) 



/o v I Mi-0 2 + /*2C 2 + /o 1 <W-C) 2 i'] 

p r is the common radial quark momentum. The parameter £ defines the position R of the 
center of mass: R = Qx\ + (1 — C,) x 2i and L is the orbital angular momentum of the system. 
The auxiliary fields fi\ and fii are seen as effective masses of the quarks whose current masses 
are m\ and m^- The last auxiliary field, v, can be interpreted in the same way as an effective 
energy for the string whose "static" energy is ar. One can get rid of these auxiliary fields 
by a variation of the Hamiltonian (J2J). Their extremal values, denoted as fiio and i/q, are the 
solutions of 

= 0, (3a) 



5H(fi h u) 
8 fii 



l^i=^i 



= 0. (3b) 

ov 

u=u 
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The equations of the RSM can be derived from the Hamiltonian (J2J by the elimination 
of v thanks to the condition (j3b|) 4j. We will consider in this paper the symmetrical case, 
where m\ = m 2 = m and /ii = \ii = fi: 



ar 2 ar Ay 2 



1- -r ^(arcsin?/ - y^l - y 2 ), (4a) 



H RS ([i) = — h fi H arcsin y + fiy 2 . (4b) 

In the general case, where m\ ^ m 2 , a third equation has to be taken into account, which 
expresses the cancellation of the total momentum in the center of mass frame p. If one 
considers fi as a number, one can directly solve Eqs. (jlj) and find the meson mass, after a 
minimization of this mass with respect to \i jj]. 

Since the RSM equations contain the remaining auxiliary field \i and a variable y whose 
physical interpretation is not a priori clear, it appears interesting to go a step further and 
to get rid of /i. As it is shown in Ref. p|, the elimination of fi with the condition Q leads 
to the extremal value 



Moreover, the replacement of /i by fio in Eqs. (Jl} leads to the following expressions 

— a = —v±-f±W r + (6a) 



tf RS M = # RFT = 2 7 ±^ r + ar aTCSmV± , (6b) 



where we have defined 



y = v ± , 1± = - / =L=, f( Vx ) = ^^^-—— and W r = y/tf + m*. (7) 

Eqs. (JHI) are precisely those of the RFTM as they appear in Ref. [| . The mysterious variable 
y is now simply interpreted as the transverse velocity Vj_ of the quarks, and the physical 
content of /i is clarified. Using definitions (J7J), we can rewrite (0) in the form 

Mo = W rlL . (8) 

Originally, the RFTM was built on phenomenological arguments. But our derivation of 
the RFTM shows that it can be derived from the Nambu- Goto Lagrangian. We also clearly 



see the equivalence of the RSM and RFTM, when the auxiliary fields are eliminated. Let 
us note that this equivalence is also true in general, with rri\ 7^ m 2 
By application of the usual correspondence rules 

1 d 2 



Pl^-- r ^r 2 r and L -> y/HJ+\), (9) 
the quantized equations of the RFTM are given by jl] 



{v±7±,W r } + a{r,f{v ± )}, (10a) 



H^ = { lL ,W r } + ^-{r,^^}. (10b) 
2 v± 

The anticommutators {A, B} = AB + BA arise because v±, r, and p r are non commuting 
operators. 

The quantized equations of the RSM, which are not used here, are given in Ref. [10]. The 



diagonalization of Hamiltonian if RFT directly provides the p 
tions (fTU|) can be numerically solved, as it is done in Refs. 



lysical masses M RFT . Equa- 
lpl. Supplemented by an 




appropriate short range potentials, like a Coulomb term, the quantized RFTM can rather 



well reproduce the meson spectra 



III. STRING AS A PERTURBATION 

When I = 0, the RFTM reduces to a spinless Salpeter Hamiltonian (SSH) with a linear 
confinement potential ar If £ is small, the contribution of the string is also small, and 
then it can be treated as a perturbation of the SSH. We will show that the contribution of 
the string can be obtained by two different procedures, that lead to different definition of 
the auxiliary field /^o- 

Let us start with the RSM and consider that the transverse velocity of the quarks is 
small: j/< 1. We can develop formulas at the second order in y. We obtain then 



H RS (n) « + n + ar + + M ) y\ (lib) 
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Inserting (jllajl in (jllbjl . and introducing p 2 = p 2 + £(£ + l)/r 2 , we can write down an 
approximate Hamiltonian 

ttA/ . p 2 + m 2 aiU + l) , in . 

defining what we call here the perturbative flux tube model (PFTM). We see that H A {[i) is 
the sum of an usual SSH with linear confinement in the auxiliary field formalism (^J 

# ss Gu) = ^±^ + /i + ar, (13) 

and a specific contribution of the string, which reads 

AH , tlM = _ jf±IL. (14) 
rp,{b/jL + ar) 



Since we want to treat Aif str (/i) as a perturbation, we will eliminate /i with the condition (|3*al 
applied for the Hamiltonian H ss (fi). This leads to the extremal value 



Ai = a/p 2 + m 2 , (15) 
and after replacement in (fT3*|) . to the spinless Salpeter Hamiltonian 



H bb = 2Vp 2 + m 2 + ar. (16) 
Let us note that formula (1151) is different from the extremal value (IS1) . The string correction 

n 

to compute is then given by [11] 

AM str = (AH st M) ~ a J+ l_KW±_ _ ; (17) 

(y/p 2 + m 2 ) ( 6(a/p 2 + m 2 ) + a(r) 



in which the mean value is performed with an eigenstate of H ss . This string contribution 
was previously obtained in Ref. ^j], where it is shown that its accuracy is better than 3%. 

We can now invert the order of the operations: Firstly to eliminate fi in the RSM and 
obtain the RFTM, then make the same approximation as before. When !)i<l, the RFTM 
equations can be developed at the second order in v± and we have 



H RFT w 2W r + ar + v 2 ± [ W r + —J . (18b) 
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Replacing Eq. (jTHaj) in Eq. ()18b|) leads to the Hamiltonian 



H A = 2W r + ar + (19) 
r 2 (W r + f) V ; 



Using the fact that 



y/P + m 2= J^ + W+ll, (20) 



we can rewrite Hamiltonian ([T9"]) in the form 



H A = 2V?W + ar + 2W r - 2\lw? + + 9 l ( f y +1) arV (21) 



/•- ' r 2 (W r + f)' 



Now, we make a new approximation and assume that W r > y/£(e + i)/r. In this case, 
which is justified in the limit v±_ <C 1, a first order expansion leads to 



^ + «£+i> M + f±i>. (22) 

r 2 2r 2 H / r 
We finally get from H A the Hamiltonian 

# A = 2Jp + mi + ar a ^+ l ) / 23 n 

rH/ r (6W r + ar) V ; 

The approximate Hamiltonian (]2Hj) can again be seen as the sum of the Hamiltonian H ss 
(jlfij) and a contribution of the flux tube 

AH - ai{i+1) (24) 

AHrit ~ ~rW r (QW r + arY (24) 

As for formula ([17]). the flux tube contribution is given by 

am, = <azw = "('+'><y> , — r , (25) 



(V^r + m2 ) (6(V^r + m2 ) + a ( r ) 



in which the mean value is again performed with an eigenstate of H ss . We immediately 
see that formulas ([T7]) and ([25]) are different. In Sec. El we will study the qualities of both 
corrections. 



IV. A MASS FORMULA 



In the following, we will specially focus on the massless case, m = 0, for which we can 
expect the largest contributions of the flux tube (/j, increases with the quark mass). In this 
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case, it is possible to find an approximate analytical mass formula. This will allow a better 
understanding of the effects of the QSE (see Sec. IVI C|) . Starting from the Hamiltonian (|12p. 
we define dimensionless conjugate variables x and q by the following scaling 

and p= (fia) 1/3 q. (26) 

We can write the Hamil- 

2\ i/3 



(/la) 1 ^ 

In this section, we treat the auxiliary field \i as a number P, ^] 
tonian (|12jl in the form 

H SS (f,)=(^y\q 2 + x)+» 
with x — \x\. Its eigenvalues are consequently 

/ a 2\ 1/3 

M ss (fi;n£) = ( — J e^ + zx 



(27) 



(28) 



where e n ^ is an eigenvalue of the Hamiltonian (q 2 + x), which is easy to solve numerically. 
Simple analytical approximate expressions of e n £ can be found for e n0 5] and e e [3] 



3tt / 3 
— 2n + - 
4 V 2 



2/3 



2 2/3 



1/3 



r(£ + 2) 
r(£ + |) 



2/3 



(29a) 



(29b) 



A more complicated approximate formula exists also in the general case Q|. 

Assuming that the quantities e n i are known, we can compute the extremal value of \x by 
a minimization of relation (|28|h which leads to 

'6n£\ 3 / 4 



A*0 nl 



3 / 



(30) 



and 



M n s f s (/i „f) = 4/i 0n £. (31) 

We will now drop the explicit dependence in n and I of the different terms to simplify the 
notations. The perturbation theory tells us that the total energy is 

ai(£+l) 



M A (/x ) = 4/i - 



(32) 



(r)// (6// + a(r)Y 

where (1/r) is replaced by l/(r). The next step is to use the Hellmann-Feynman theorem 
E Q which states that 

M ss (/i ) =4/i = 2/i + a(r). (33) 
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Extracting (r) from relation (j33J) and replacing it in relation (|32jl . we finally obtain the mass 
formula 

M^)^A^- aH [t\ l \ (34) 

where /xo is given by formula (jHOJ)- 

The mass formula ()34j1 . even approximate, exhibits Regge trajectories. This can be easily 
checked when n = 0. Using the fact that 

lim e , « 3 — ± , (35) 



2 

with formula (|30j). we obtain the following expression for the extremal value of the auxiliary 
field at large values of £ 

fa£ 

/*>«Vy (36) 

Replacing Eq. (jHEJ) into Eq. pijl . we get 

(M A M) 2 -f^. (37) 

This linear relation between the squared mass and the angular momentum reproduces qual- 
itatively the Regge trajectories. The Regge slope is here 7.03 a, a higher va,ue than the one 
predicted by the RFTM, which gives a slope equal to 27ia |1| . Finally, we can observe that 

lim = -. (38) 

e^oo 4/i 16 

When I = 0, this ratio is vanishing. This justifies to treat the contribution of the string as 
a perturbation. 



V. COMPARISON BETWEEN EXACT AND APPROXIMATE STRING CON- 
TRIBUTION 

Before studying the contribution of the QSE to the RFTM, it is interesting to examine 
the relevance of the approximate treatment for the flux tube developed in Sec. IIHI For this 
purpose, we compare here some masses computed with the "exact" RFTM by numerically 
solving Eqs. (jlOj) . as it is done for example in Ref. ^(|, with masses computed in the 
framework of the PFTM with both string corrections (|17J) and (|25|). 
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The symbol Me designs a mass computed with the RFTM, and Mpi, Mp 2 the corre- 
sponding PFTM masses evaluated with the string corrections (ITT)) and (J25J) respectively. 
The quantities 



AM, 



E,Pi 



M E - M Pi 



M f 



(39) 



measure the differences between the RFTM and its approximations. 



TABLE I: AMe,pi (correction (|17|)) in %, for different states, with m = 0. £ is the orbital angular 
momentum of the state and n is the number of nodes at finite distance. 



n = 
1 
2 
3 
4 
5 



1 

3.26 
1.08 
0.68 
0.46 
0.36 
0.29 



3.78 
1.31 
0.86 
0.56 
0.43 
0.35 



4.01 
1.53 
1.04 
0.68 
0.52 
0.40 



4.17 
1.74 
1.22 
0.81 
0.64 
0.49 



4.29 
1.94 
1.40 
0.95 
0.76 
0.58 



Table |l] presents the difference between the RFTM and the PFTM with correction (|17|). 
As the basis of the approximation was to consider a small v±, it is not surprising to see 
that AMe,pi increases with £. On the contrary, AMe,pi decreases for an increasing radial 
quantum number n. This can be understood by the presence of the operator p% in the 
denominator of Eq. ()17j) . This term becomes larger with n, and makes the contribution 
of the string smaller, leading to a decreasing of AMe pi- Globally, the approximation 
considered in Table H] is rather good, in particular when n > i. 

Table ITT1 presents the difference between the RFTM and the PFTM with correction ()25|) . 
When n = 0, we immediately see that this approximation is not so good, because AMe,p2 
becomes large very quickly with £. The situation is better for larger values of n, when 
our approximation W r ^> \J l{l + is particularly justified. But the evolution is less 
monotonic than in Table UJ In particular, for a fixed value of n, AMe,p2 decreases to a 
minimal value, and then increases again with £. In conclusion, in the framework of the 
PFTM, the correction (fTTj) . proposed in Ref. , seems preferable (even if the other one 
can sometimes give a better result). It provides quite globally good results. 
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TABLE II: Same as in Table [Qbut for AM EtP2 (correction (12310 . 



n = 1.03 7.81 15.21 22.52 29.54 

1 0.44 0.88 2.55 4.41 6.36 

2 0.40 0.67 0.73 1.62 2.64 

3 0.35 0.14 0.23 0.72 1.30 

4 0.31 0.19 0.02 0.30 0.66 

5 0.26 0.19 0.06 0.13 0.37 

VI. THE QUARK SELF-ENERGY 
A. Definition 

Recently, it was shown that the QSE contribution, which is created by the color magnetic 
moment of the quark propagating through the vacuum background field, adds a negative 
constant to the hadron masses 8fl. Its negative sign is due to the paramagnetic nature of 
the particular mechanism at work in this case. Using the Fock-Feynman-Schwinger repre- 
sentation of the quark Green's function, one can obtain the QSE contribution as a shift of 
the squared mass of the quark P, [3] which reads 

POO 

Am 2 = -3m / dzz 2 K x {mz) (D(z) + D^z)) , (40) 
Jo 

where D and D\ are quark correlators and K\ the Mac-Donald function. The properties 
of these correlators were studied by lattice simulations in the quenched case One has 
then 

D(z) xsSD^z) = D(0)exp(-\z\5), (41) 

with 5 = 1/Tg. T g is the gluonic correlation length, whose value is estimated at about 
0.15-0.2 fm. This locates 5 in the interval 1.0-1.3 GeV. The results (j41j) allow us to find an 
analytic form for the integral (jjOj) 

Am 2 = -4mD(0)cp(m/5), (42) 
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where, defining e = m/8, we can write ip(e) as 



1 


-3e 


0° 


(1 - e 2 ) 7 


1 


-3e 


53 


(e 2 _ l} 5/2 



mi 1 + ^V ^ 



arctan (\/e 2 — l) 



e(l -e 2 ) 2 
l + 2e 2 



<5 3 



(<0 (c < 1) 



e 1-e 



2^ 



(43) 



5 3 



One can check that 0i(l) = 02(1) = 2/5, and that 

lim0i(e) 



(44) 



For a purely exponential correlator, as it is the case here, D(0) is connected with the 
string tension a by the relation 



so we find 



2 4am5 
Am = V 9 ( e )- 

7T 



(45) 



(46) 



For convenience, we define a new dimensionless function, 77(e), by 

77(e) = 5 3 e^(e). 



(47) 



We see in Fig. [J that 7/(0) = 1 and that 77(e) rapidly decreases for increasing values of e. 

1.0 




FIG. 1: Plot of 77(e). 

Thanks to the definition (|47J). Eq. ({US]) takes its final form 

Am 2 = 77(e). 

7T 



(48) 
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An important ingredient we have used to get the contribution (|48|) of the QSE is the rela- 
tion (jUJ, derived in the quenched case. In Ref. the results obtained in the unquenched 
case are quite different: the exponential form of D remains the same, but now D\ is small 
enough to be neglected. This approximation leads us, after the same calculations as before, 
to 

Am 2 = -— 77(e). (49) 

7T 

The two formulas (|48p and (J4U|) only differ by a constant factor: 4 in the quenched case, 3 
in the unquenched one. Since this factor does not seem to be presently known with a great 
accuracy, we will finally use the following expression for the quark self energy 

Am 2 = -—rj(m/5), (50) 

7T 

with / e [3,4] and 5 G [1.0, 1.3] GeV. 



B. Insertion of QSE in effective meson models 

In the previous section, we showed how the QSE contribution acts as a shift of the squared 
mass of the quarks. We have now to insert this new term in the models we described in 
Sees. ITT1 and ITTT1 If we make the substitution m 2 — > m 2 + Am- in the Hamiltonian (j2J, we 
find 

H - H + AH QSE , (51) 

where 

a tt Am i f a V( m i/ $) /, 9 x 

Equation ()52j) is the total contribution of the QSE to the RS Hamiltonian. This term has 
to be considered as a perturbation of the original Hamiltonian, and thus one has not to give 

n 

much sense to the fact that for light quarks the total mass m 2 + Am 2 is negative |SJ. Since 
the QSE is a perturbation of the Hamiltonian, it has not to be included in the elimination of 
the auxiliary field \l. To take into account the QSE in the RFTM, we suggest the following 
procedure, inspired from Ref. 8j: 

1. To find the eigenvalues and eigenf unctions of the RFTM Hamiltonian. 

2. To compute the mean value (/j,q) of the extremal field ^0 with the eigenf unctions. 
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3. To add to each eigenvalue the corresponding QSE contribution which reads, in 
the symmetrical case, 

far)(m/S) 

AMqse = -— A-T 2 - 53 

The problem is to choose the value /io of the extremal field. Within the PFT, this value 
is given by Eq. (|15)) . and the resulting QSE correction is given by 



7T (a/p 2 + m 2 ) 

as it is done in Ref . . On the other hand, if the equations of the RFTM are not treated 
in perturbation, it seems natural to take W r , y± for /i (see Sec. |H}. Within this framework, 
the QSE contribution is expected to be given by 



A7lf RFT fa y(m/5) 

AM QSE = (55) 



Actually, we have to replace W r j± by {W r ,7j_} /2 in order to keep the operator hermitian. 
We can expect that both procedures will lead to different results, as in the case of the 
contribution of the string as a perturbation (see Sec. IHIjl . 

C. Regge trajectories 



As mentioned in Ref. p|, the QSE correction preserves the Regge trajectories. We can 
qualitatively understand this thanks to the mass formula (|34|). to which we add the QSE 
contribution (J52J) . We have 



mv,)= ^_^^_fjL^m, (56) 

16/ijj 7T /i 



For large angular momentum, /x becomes large. Keeping only the dominant terms, we find 
the approximate mass formula 

(M A (/i )) 2 ^16^-8^(m/5). (57) 

7T 

It appears that the Regge trajectories are preserved, since the QSE only causes a global 
shift of the squared masses and preserves the dominant 16/Xq term, which grows like i. This 
is particularly clear when n = (see Eq. (J3BJ)). 
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VII. ADDING THE QSE 



In this section, we compute numerically the contributions of the QSE for both the RFTM 
and the PFTM. We will focus on the massless case, for which we expect the largest effect 
since rj(0) = 1 and the extremal field fio is minimal. In this special case, we are able to 
perform a universal analysis of our results, since the meson mass are then just scaled by 
the factor y/a. We will use in this section / = 3.0, which is the value computed in the 
unquenched case. 




FIG. 2: Main : Regge trajectories for n = with (circle) and without (triangle) QSE, computed 
with the RFTM. Small box: Same Regge trajectories computed with the PFTM. Lines are used 
to guide the eyes. 

We noticed in Sec. IVI CI that adding the QSE contribution did not destroy the Regge 
trajectories and the Regge slope. This result was obtained using a mass formula, itself an 
approximation of the PFTM. Since we are able to numerically solve the RFTM without 
making approximations Q], we can directly study the influence of the QSE correction on 
the Regge trajectories. The main graph of Fig. El immediately shows the negative shift 
of the squared masses when the QSE is added, as expected. If the linearity of the Regge 
trajectories is well preserved with the RFTM, the Regge slope is not. In this figure, j3 is 
the Regge slope with QSE and a is the corresponding one without QSE. Both are rather 
different, and we obtain (3 /at — 0.88. This diminution of the slope is caused by the QSE 
term (|55p. The small box shows that when one is working with the PFTM and the QSE 



15 



term (|54*j). as it is the case in [llj], the Regge slope is not affected, or a very little bit. We 
find indeed a = 7.00 and (3 = 6.95. 



m= 1.0 GeV 
m = 0.6 GeV 



4 6 
/ 



FIG. 3: QSE contributions, for n = 0, versus t. The QSE contributions in the RFTM (filled circle) 
and with the PFTM (open circle) are computed for different quark masses. Note that when m ^ 0, 
the scaled results are not independent of a, which is taken here equal to 0.19 GeV 2 . 

We have now to understand why the QSE affects the Regge slope in the RFTM and not 
in the PFTM. Fig. El illustrates the differences between the QSE contributions ()54)1 and (|55jl. 
When i = both contributions are equal since the PFTM and the RFTM reduce to the same 
spinless Salpeter equation. In the massless case, the two contributions considerably differ 
for £ ^ 0. The QSE coming from the RFTM is always the smaller one, and this causes the 
Regge slope to be smaller when one solves the RFTM. As expected, the difference between 
the two contributions decreases when the quark mass increases. With a large quark mass, 
both RFTM and PFTM have a common non relativistic limit (when the approximation 
v j_ <C 1 is the most justified). The conclusion to draw from Fig. El is that (W r r y± ) is not 
even approximately equal to ( ^/p 2 + m 2 } for light quarks. Moreover, the second expression 
leads to a QSE term which preserves the Regge slope, while the first one does not. We show 
in Fig. |U the effect of adding to the RFTM solutions a "theoretically justified" QSE term 
(155(1 and a "PFTMdike" one ([54)1 . As expected, the contribution (p>4*|) causes no diminution 
of the Regge slope. 
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FIG. 4: Regge trajectories, with n = 0, for the genuine RFTM (filled circle), the RFTM with 
a theoretically expected QSE term AMq^g (open circle) and with a QSE term inspired by the 
PFTM AM<^ (triangle). 

VIII. EXPERIMENTAL DATA 



Since the RFTM includes neither the spin (S) nor the isospin (J) of the mesons, the 
experimental data we will try to reproduce here are the spin and isospin averaged masses, 
denoted M av . These are given by [17] 

E,„,(2I+1)(2J + 1)M„ 
" E,.,(2/ + D(2J + 1) ' (58) 

with J = L + S and M/ j are different masses of the states with the same orbital angular 
momentum I. The first three columns of Table IIVI show the experimental data concerning 
nn {n means u or d) states used to compute the spin-isospin averaged masses. These data 
are taken from Ref. jis| . 

Realistic masses cannot be computed with the genuine flux tube Hamiltonian (jlObjl . It 
can indeed give the right Regge slope with a < 0.2 GeV, but the absolute values of the 
masses are always too high. An attractive Coulomb potential, simulating the one-gluon 
exchange process, must be added 

V(r) = (59) 

where as is the strong coupling constant. Theoretical arguments as well as lattice calcula- 
tions agree with a value as < 0.4 19]. However, even such a Coulomb term does not shift 
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the masses enough to reproduce the data (see for example Ref. jjjj), except if unrealistic too 
high values are chosen for as- This means that other contributions coming from neglected 
physical mechanisms are still needed. The QSE seems to be an interesting one, because its 
contribution to the mass is negative, and rather large. This should allow to reproduce the 
experimental data with an acceptable value for the Coulomb term. 

TABLE III: Two sets of possible physical parameters. 



Type 1 Type 2 

a (GeV 2 ) 0.245 0.195 

m n (GeV) 0.0 0.073 

a s 0.4 0.4 

/ 3.0 3.0 

S (GeV) 1.0 1.0 

QSE term AM^ AM^ 



In Table IHIl we give the two sets of parameters we use to compute the masses of some nn 
states. Both have the same value of /, 5 and as but differ for the other parameters. In the 
type 1, we make the usual choice m n = and take formula (|HH|) as QSE term. This is the 
contribution that one can theoretically expect. As we have seen that it causes a diminution 
of the Regge slope, we have to take a = 0.245 GeV 2 in order to obtain a final slope in 
agreement with the experiment. Choosing a = 0.2 GeV 2 is no longer possible as it is the 



case with the PFTM 



11 1 . This is the unconventional aspect of type 1 set of parameters. 



On the other side, with the type 2 set, it is possible to keep for a the standard value, about 
0.2 GeV 2 . The price to pay is to take formula (|54jl as QSE term, an only "empirically" 
justified choice, and to give a small mass to the quark n. The comparison between the 
experimental averaged masses and our results is given in Table IIVI We see that both types 
lead to masses close to the spin-isospin averaged ones, our results being located inside the 
error bar in almost every case. 
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TABLE IV: Comparison between the spin averaged masses M av of the nn states and the results 
of the RFTM plus a QSE term. M\ and M 2 are the masses computed with the set of parameters 



1 and 2 respectively. 


Masses are given in GeV. The first three columns show the different states 


used to compute 


the 


spin averaged masses. 




State 


I 


(n + l) 2S+1 Lj 


M av Mi M 2 







i 3 s 1 


0.773 ±0.011 0.788 0.772 


P 


1 


I 3 Si 




fci(1170) 





I'Pi 


1.265 ±0.011 1.269 1.282 


6i(1235) 


1 


I'Pi 




/i(1285) 





1 3 P 1 




01 (1260) 


1 


l 3 P 1 




/ 2 (1270) 





1 3 P 2 




a 2 (1320) 


1 


1 3 P 2 




w(1650) 





1 3 D 1 


1.676 ±0.012 1.673 1.678 


p(1700) 


1 


l 3 D x 




w 3 (1670) 





1 3 D 3 




p 3 (1690) 


1 


1 3 D 3 




/ 4 (2050) 





1 3 F 4 


2.015 ±0.012 2.016 2.006 


a 4 (2040) 


1 


1 3 F 4 





IX. CONCLUDING REMARKS 



The purpose of this work was to study the contribution of the quark self energy to the 
meson masses in the framework of the relativistic flux tube model. The equivalence between 
this phenomenological model and the more QCD based rotating string Hamiltonian is used 
as a guide to perform the calculations. 

The equations defining the relativistic flux tube model being rather complicated to solve, 
it seems interesting to treat the flux tube contribution as a perturbation. Two procedures 
have been studied: To eliminate firstly the auxiliary field from the rotating string Hamil- 
tonian, or to make firstly the approximation of small transverse velocities in the rotating 
string Hamiltonian. We arrived in Sec. IIHI at two non equivalent terms for the flux tube 
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correction: the first one, obtained in the auxiliary field formalism, was already known [11], 
and the second, obtained directly from the relativistic flux tube equations, is a new one. 
The results of both contributions are quite similar, but we showed in Sec. E] that the first 
approach is globally the best one. 

Starting from the rotating string Hamiltonian and considering the auxiliary field asso- 
ciated with the quark mass /i as a simple number, an approximate but analytical mass 
formula is established. It enables to understand at least qualitatively why the quark self 
energy preserves the linearity of the Regge trajectories and decreases the squared masses by 
a constant quantity (see Sees. IIVI and IVI (If . 

The addition of the quark self energy to the relativistic flux tube model (Sec. IVI 1)1 pre- 
serves the linearity of the Regge trajectories. But, for massless quarks, the Regge slope is 
smaller by a factor 0.88 with the quark self energy than without it (this value tends toward 
unity when the quark masses increase) . This effect does not exist when one works within the 
framework of a perturbation theory, in which the relativistic flux tube Hamiltonian reduces 
to a spinless Salpeter Hamiltonian. It is due to the fact that different extremal values of the 
field fi are found, y/p 2 + m 2 or W r 7±, considering the perturbation scheme or not. 

In the framework of the relativistic flux tube model, it is possible to reproduce the 
experimental data with the theoretically expected quark self energy term (// = W r , y±) and 
an realistic Coulomb term, if larger value than usual, 0.245 GeV 2 , is chosen for the string 
tension a. It is possible to keep the well known value a = 0.2 GeV 2 by using a quark self 
energy term coming from a perturbation approach (// = \/p 2 + m 2 ). This leads us to two 
possible conclusions. Within the relativistic flux tube model: 

• It is necessary to use a quark self energy term in which the extremal field \x is given 
by y/p 2 + m 2 in order to keep a value of the string tension around 0.2 GeV 2 . 

• It is necessary to use a quark self energy term in which the extremal field /i is given by 
W r r y±, the natural value associated with the relativistic flux tube model; But another 
physical mechanism has to be taken into account in order to keep a value of the string 
tension around 0.2 GeV 2 . 

It seems hard to justify the first proposal by any theoretical argument. The mentioned new 
mechanism could be due to retardation effects or deviations from the straight line for the 
flux tube. Such a work is in progress. 
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